Optimal Weighting in Galaxy Surveys: 
Application to Redshift-Space Distortions 
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Using multiple tracers of large-scale structure allows to evade the limitations imposed by sampling 
variance for some parameters of interest in cosmology. We demonstrate the optimal way of carrying 
out a multitracer analysis in a galaxy redshift survey by considering the principal components of 
the shot noise matrix from two-point clustering statistics. We show how to construct two tracers 
that maximize the benefits of sampling variance and shot noise cancellation using optimal weights. 
On the basis of high-resolution TV-body simulations of dark matter halos we apply this technique 
to the analysis of redshift-space distortions and demonstrate how constraints on the growth rate of 
structure formation can be substantially improved. The primary limitation are nonlinear effects, 
which cause significant biases in the method already at scales of < 0.1/iMpc~^, suggesting the 
need to develop nonlinear models of redshift-space distortions in order to extract the maximum 
information from future redshift surveys. 

PACS numbers: 98.80.-k, 98.62.-g, 98.65.-r 



I. INTRODUCTION 

One of the deepest mysteries of contemporary cosmol- 
ogy is the nature of the observed accelerated expansion 
of the universe. So far, its evolution can be described 
remarkably well by Einstein's theory of gravitation in- 
cluding a nonzero cosmological constant A. However, in 
order to fit the astronomical observations (e.g., pj]), A 
must be many orders of magnitude smaller than what our 
standard model of particle physics would expect. This 
hierarchy-problem inspired various departures from the 
cosmological standard model, such as modifications of 
Einstein's field equations or the introduction of exotic 
forms of matter. 

A particularly sensitive probe of cosmology is the large- 
scale structure (LSS) of the universe. Galaxy redshift 
surveys map out large fractions of its observable vol- 
ume and thereby reconstruct a three-dimensional map 
of density fluctuations, whose statistical properties di- 
rectly relate to fundamental cosmological parameters [ll- 
Unfortunately, this reconstruction is hampered by 
the fact that galaxies are biased and stochastic tracers 
of the dominating dark matter density field. Even if the 
density field could be inferred perfectly well, the finite 
number Nk of independent Fourier modes in the survey 
sets a fundamental lower limit on the achievable uncer- 
tainty, which is known as sampling variance (or cosmic 



* IhamausQphy sik. uzh . ch | 



variance, in case the survey size is the whole observable 
universe). For example, in a measurement of the dark 
matter power spect rum P , the sampling variance limit is 
given by ap/P > \plJWk, an uncertainty floor that prop- 
agates into all the parameters one wants to infer from P. 
This limit decreases towards smaller scales as more and 
more Fourier modes can be sampled, but, at the same 
time linear theory starts to break down and higher-order 
perturbation theory has to be adopted to model P (see, 
e.g., 0-11] )• Further complication arises in relating the 
observed galaxy power spectrum to the latter, as galaxy 
bias becomes nonlinear and nonlocal [13] . 

An alternative approach to accurately probe cosmo- 
logical parameters is to consider multiple tracers of the 
density field within the well-understood linear regime. 
The relative clustering amplitude between multiple trac- 
ers can be inferred without sampling variance limitation, 
because the underlying density fluctuations cancel out in 
taking ratios 12]. Therefore, any cosmological infor- 
mation that remains in the relative clustering amplitude 
between different tracers can potentially be inferred with 
a much higher accuracy. 

In this paper we focus on a particular contribution to 
the clustering amplitude of galaxies coming from redshift- 
space distortions (RSD). These are caused by peculiar 
velocities along the line of sight, causing their clustering 
statistics to become anisotropic. First treated as a con- 
tamination, this effect has been realized to be a powerful 
probe of cosmology, as an understanding of RSD allows 
to infer the growth rate of structure formation, which is 
directly tied to the expansion history of the universe as 
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well as the theory of gravity (e.g., [13413 )• 

After recapping the fundamentals of RSD and intro- 
ducing a general formalism for the multitracer analysis 
in Sec. HIl we present our results on the RSD analysis 
from TV-body simulations in Sec. IIIIl Finally we draw 
our conclusions in Sec. IIVI 

II. FORMALISM 

A. Galaxies in redshift-space 

In galaxy surveys, radial distances are inferred via the 
individual redshift of objects, assuming they follow the 
Hubble flow. Due to gravitational attraction however, 
galaxies (respectively their host halos) build up pecu- 
liar velocities v which contribute to their redshift via the 
Doppler effect. Hence, the real-space and redshift-space 
locations r and s of a galaxy are related as 



where r is the unit vector along the line of sight and H{z) 
is the Hubble constant as a function of redshift z. On 
large scales, gravity causes coherent infall of test particles 
into the potential wells of the dark matter. Thus, galax- 
ies are moving towards overdense regions in the universe, 
resulting in an enhancement of their inferred overdensity 
along the line of sight. According to linear perturbation 
theory, the redshift-space and real-space galaxy overden- 
sities are related as 

<5(^)(k,M)=4'-)(k) + /Ai2<5(k), (2) 

where /i the cosine of the angle between any wave vector 
k in the survey and the line of sight, 6 the dark matter 
overdensity field and / the growth rate of structure. This 
well known result (2^] further makes use of the plane- 
parallel approximation, assuming the separation of any 
galaxy pair to be much smaller than their distance to 
the observer. On nonlinear scales, random motions are 
generated in the process of virialization, which causes a 
damping in the clustering amplitude of galaxies along 
the line of sight. This so-called Finger-of-God (FoG) 
effect is often modeled phenomenologically by an addi- 
tional Gaussian damping factor in Eq. ^ [25|] , but more 
elaborate schemes have been developed (e.g., p6l - l33 |). 
We will neglect nonlinear corrections in this paper and 
focus on the large linear scales. 

The growth rate / is the logarithmic derivative of the 
growth factor D with respect to the scale factor a. In 
linear theory, it can be expressed as 

/ = dlnZJ/dlna ~ 172, , (3) 

with the matter density parameter and the growth 
index 7 [ssj . In Einstein gravity, for most models of 



dark energy the value of 7 is about 0.55, but can take 
on distinctly different values in modified gravity scenar- 
ios (3^. Therefore, constraints on the growth rate can 
provide viable tests on the theory of gravitation. 

Galaxies only form in specific, discrete locations of the 
density field, they are referred to as biased and stochastic 
tracers of the dark matter. On linear scales, this relation 
can be described locally by 

S^;\r)^b,S{r) + e, (4) 

where the factor bg is the linear galaxy bias and e a 
random variable denoted as shot noise, describing the 
stochastic nature of this relation. Both bg and e depend 
on redshift, as well as various properties of the type of 
galaxies one is considering (e.g. luminosity, color, host- 
halo mass). 

Together with Eq. ([2]) this yields a model for the over- 
density field of galaxies in redshift-space, 

4^)(k,Ai) = (6g + V)5(k)+6. (5) 

Since the phenomena of galaxy biasing and RSD are to 
multiply the density field 6 by some factor, we can sim- 
ply define a more general effective bias parameter b that 
contains both contributions, 

b = bg + f^i^ . (6) 

In the following, we will drop the superscript that dis- 
tinguishes between real- and redshift-space quantities for 
clarity. If not explicitly mentioned otherwise, all symbols 
should be understood as given in redshift-space. 

B. Multiple tracers 

In order to exploit the gains of sampling variance and 
shot noise cancellation, we need to consider multiple trac- 
ers of the dark matter density field [ll|j [13 • One way to 
achieve this is splitting some galaxy catalog into bins of a 
certain observable property of the galaxies (like luminos- 
ity, color, host-halo mass, etc.). However, the following 
framework is not limited to galaxies and may be adopted 
for other tracers of the dark matter density field as well. 

1. Covariance matrix 

We start by writing the density fields of N tracers as a 
vector 5g = ((5g, , (5g2 , . . . , ) . The outer product of this 
vector, once ensemble averaged within a fc-shell in Fourier 
space, yields the covariance matrix C = (5g(5g), where 
the ^^-symbol denotes the hermitian conjugate (transpose 
and complex conjugation). Plugging in the model for 
galaxy overdensities in redshift-space from Eq. ([5]), it 
reads 

C = bb''P + £, (7) 
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with the effective bias vector b = bg + //i^, the dark 
matter power spectrum P = {66*) and the shot noise 
matrix £ = (ee^) (by definition {e6*) — 0). The covari- 
ance matrix contains all auto- and cross power spectra 
of the considered tracers, so in total N{N — l)/2 inde- 
pendent elements per /c-shell. However, since all tracers 
follow the same dark matter density distribution, these 
elements are correlated. 

Let us compare the analysis for a single tracer with 
the one for two tracers. From a single tracer we can only 
observe an estimator of its auto-power spectrum 



(8) 



where we define /3 = f/bg. In this parametrization it 
is obvious to see that bg is degenerate with the power 
spectrum P, respectively its normalization as defined via 



(9) 



where Pq describes the shape of the power spectrum and 
is assumed to be known, at least up to linear order. From 
Eq. ([5]), only the combination b^P can be determined at 
/X = (usually, the shot noise is assumed to be Poisso- 
nian, meaning it is scale-independent and given by the 
inverse number density of galaxies n~^). The same is 
true for the growth rate /, we can only determine the 
product 



C{k, n)-£- [C{k, 0) -£]{! + 2/3/^2^ 



(10) 



cancel out completely and the three ratios will be dif- 
ferent and scale-dependent. This gives rise to residual 
sampling variance inherent to a measurement of a and 
/3, which in general can be much smaller than in the 
single-tracer case. In turn, a better measurement of /3 
allows a more precise estimate on /^P (see Eq. ([T0|'). 
However, the accuracy on f^P is still limited by sam- 
pling variance, yielding (Tf2p/pP > ^ 2 jN-k , respec- 
tively a^^Jjos > ^/T/2N; 

If both and 6 are known, we can simply add the 
dark matter overdensity mode S to the overdensities of 
the tracers and write S = {6,6g). In this case, the ef- 
fective bias can be obtained directly by taking the ratio 
{Sg6*) / (SS*) . The covariance matrix then becomes 



C = 




P b^P 
bP bb^P + £ 



(15) 



with N{N + l)/2 independent elements. Now the degen- 
eracy between &g, erg and / is lifted, because P is known 
separately. Considering the cross-correlation coefficient 
of a tracer Sg with the dark matter 6, we find 



{6g6* 



{Sg6*g){66* 



= 1 



£ 



{bg + f^i^yp, 



(16) 



so the deviation of r from unity is only controlled by the 
shot noise of the tracer. Again, if f = 0, the dark matter 
power spectrum disappears completely, so that bg and / 
can be determined without sampling variance. 



but this does not apply for /3, which can be extracted 
directly from observations of C where fi ^ (e.g., [131). 
Note, however, that in this case the achievable error on 
f3 is limited by the sampling variance inherent to P. 

In case two distinct tracers of the density field with 
biases bg and abg are observed, where a is their relative 
galaxy bias, we obtain the three following power spectra: 

Cn{k,pi) = {l + l3^i^)HlP{k)+£n , (11) 

C22(fc,A*) - {a + l3n^)hlP{k)+ £22 , (12) 
Ci2(fc, m) = (1 + /3/^') (a + blP{k) + £12 ■ (13) 

The degree to how well they are correlated is quantified 
by the cross-correlation coefficient = C12/ CiiC22, so 
in the idealistic case of no shot noise (5 = 0), r = 1 and 
the ratios 

C12/C11 = C22/C12 ~ \/ C22/C11 — ^ ^ (14) 

all yield the same expression, which is independent of 
P. Hence, a combination of observations at different 
fj, yields a and ^ without sampling variance P]. Un- 
fortunately, in realistic surveys £ 7^ 0, so P will not 



2. Fisher information 

In order to determine more quantitatively how much 
information on cosmology is buried in the clustering 
statistics of biased tracers and the dark matter, we have 
to compute the Fisher information matrix [38], a deriva- 
tion of which is presented in the following: we start with 
a multivariate Gaussian likelihood of the data- vector ^„ 



1 



(27^)^/2^/d^^C 



exp 



1 



(17) 



which is a reasonable assumption on large scales, where 
(5 <C 1. The Fisher information matrix for the parame- 
ters 9i and 9j is obtained by ensemble averaging over the 
Hessian of the log-likelihood [li,!!^!, 
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According to the model from Eq. (O, the derivative of the 
halo covariance matrix with respect to the parameters is 



1^ = {bbl + b,b^^P + bb^P, 



(19) 
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where bi = db/d6i, Pi = dP/dOi and we assume 
d£ /dOj — 0. Utilizing the Sherman- Morrison formula 
El, , the inverse of the covariance matrix becomes 



l + b^£-^bP 



The full Fisher matrix for this case is calculated in Ap- 
pendix [X] The result is 



F — 



where 



E P 



E P 



PP, 



2P2 



E,, ^ b\£ 



-%P 



(21) 



(22) 



and Si as well as E are defined accordingly by simply 
omitting the corresponding indices (derivatives) . We can 
identify the first two terms in the square brackets as a 
single tracer- and a multitracer term, respectively. The 
single-tracer term is suppressed by a factor as com- 
pared to the multitracer term. By definition, the lat- 
ter vanishes for the case of only one tracer, since then 



£ 
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EEi, = EiEj and Eq. (|2T|) simplifies to 
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b P 



P P 
6 P ^ 2P2 



(23) 

While the Fisher information for multiple tracers can in 
principle become infinite in the limit of no shot noise, 
Eq. ([23l) reaches a finite limit when f — 0. 

We consider two sets of parameters separately: the 
ones that influence only the effective bias. 



(24) 



and the ones that go into the matter power spectrum, 

= ((78, n„ h, Qa, an, f^b, f^k) . (25) 

The elements of the shot noise matrix are usually not 
considered as quantities of interest, but as nuisance pa- 
rameters that can be marginalized over. From Eq. (|2ip 
it is evident that there are degeneracies between and 
0^^'' due to mixed terms of the form EiP,-. However, 
if E is sufficiently large, those terms are suppressed in 
comparison to the term E^j. 

In case the dark matter density field is known in ad- 
dition to the galaxies, the derivative of C from Eq. (|15l) 
becomes 



dC 



P 

biP + bp 



b]P + b''P 
tj 



bb\P + b,b^P + bb^P, 



(26) 



Furthermore, a block-inversion of C yields 



(1 



|-E)F-i -b^£ 
£^b £-^ 



(27) 



(20) and the Fisher information becomes (see appendix El) 



2P2 



For a single tracer this further simplifies to 



_ hb,P 



P.P, 

2P2 



(28) 



(29) 



This expression may increase indefinitely for sufficiently 
small £. In the limit 5 — > 0, the effective bias can be 
determined exactly, allowing an exact measurement of 
the parameters [11]. On the other hand, constraints 
on the parameters 0*-^'' are always limited by the variance 
of the power spectrum Var(P) = 2p2. Here, there are 
no mixed terms depending both on the effective bias and 
the power spectrum, so parameter degeneracies between 
e^^'> and e^^'^ are absent. We note that in the limit of 
small E.,/E, Eq. ^ reduces to Eq. p8)) . 



3. Reparametrization 

As mentioned above, the growth rate / cannot be de- 
termined from a galaxy redshift survey alone, since it is 
degenerate with the power spectrum P, respectively its 
normalization erg. This degeneracy can only be broken 
with knowledge of the dark matter density field, or other 
prior constraints. For this reason it is sometimes conve- 
nient to reparametrize by the mapping 



0''^=0(76f =(bg,/)/6f ^(a,/3) , 



as = (78 6[ 



rcf 



(30) 
(31) 



where 6g°^ is an arbitrary reference galaxy bias so that 
a is the relative galaxy bias of all considered tracers to 
this reference. This mapping leaves the covariance matrix 
from Eq. ([7]) unchanged. We can conveniently choose the 
lowest bias of all tracers as the reference, 6g°* — 5g^ , such 
that a = (1, 6g2/6gj , . . . , 6g„/6gJ is always larger than 
unity and f3 = f /b^-^ . Another popular parametrization is 



(32) 
(33) 



with (7™^ being some arbitrary reference normalization 
of the power spectrum, which we choose to be identical 
with our simulation input value of ag = 0.81. In this 
paper we will quote constraints on both /3 and /cts when 
considering a galaxy redshift survey only, and on / when 
the latter is combined with dark matter observations. 
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III. ANALYSIS 
A. Numerical setup 

Our numerical analysis is based on a high-resolution 
A^-body simulation performed at the University of Ziirich 
supercomputer ZBOx3 with the gadget ii code [43]. It 
contains ISBG"^ particles of mass 4.7 x IO^^/i^^Mq in a 
box of 1.3/i^^Gpc a side. We chose our fiducial cos- 
mology to match the wmap5 best fit with ag = 0.81, 
ris = 0.96Ji = 0.7, riA = 0.721, fl^ = 0.279, fl^ = 0.046, 
r^k = l44i . We further employ a friends-of-friends al- 
gorithm [4^ with a linking length of 20% of the mean 
inter-particle distance to generate halo catalogs at dif- 
ferent redshifts. With a minimum of 20 dark mat- 
ter particles per halo we resolve halo masses down to 
Mmin ~ 9.4 X 10"/i-1Mq. 

In order to transform the real-space halo catalog into 
redshift-space, we apply Eq. ([ij using the velocities 
of the halos along the three independent directions of 
the box (x-, y- and z-axis). Thus, three independent 
redshift-space catalogs can be constructed from a single 
real-space catalog, yielding a total effective volume of 
Ves — 6.6/i~'^Gpc'^. Density fields are created by cloud- 
in-cell interpolation 0| onto a mesh of 1024^^ grid points, 
and the Fourier-modes are obtained using an FFT algo- 
rithm. 

We utilize the idl algorithm mpfit 47] to fit our mod- 
els to the numerical data and find the best-fit parameters 
including their uncertainties. It is based on the minpack 
distribution by and uses the Levenberg-Marquardt 
technique to find the minimum of a multi-dimensional 
nonlinear least-squares problem. Parameter uncertain- 
ties are calculated via the Jacobian of the chi-square, 
which is determined numerically using finite-difference 
derivatives. 



B. Optimal tracer selection 

In order to fully exploit the benefits of the multitracer 
approach, the question remains on how to ideally con- 
struct different tracers from a given galaxy catalog. In 
this section we will derive the answer to that question 
and test it on the basis of A^-body simulations. 



1. Principal components of the signal-to-noise ratio 

A quantity that plays a crucial role for the cosmologi- 
cal information content contained in the two-point clus- 
tering statistics of LSS is the clustering signal-to-noise 
ratio S = b^S^^bP. For dark matter halos, it has been 
demonstrated on the basis of numerical A^-body simula- 
tions that E is dominated by only two principal compo- 
nents, corresponding to the two non-trivial (non-Poisson) 
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FIG. 1. Three different eigenvectors of the shot noise matrix 
obtained from a halo catalog that has been split into 30 mass 
bins of equal number density. The two non-Poisson eigenvec- 
tors (red and blue stars) are overplotted with their best-fit 
weighting functions w± as solid lines with functional form 
given in the top left of the panel in corresponding colors. One 
representative Poisson eigenvector is shown as green circles, 
when used as weighting function it yields a very low clustering 
signal-to- noise ratio S, due to its oscillatory behavior. 



eigenvectors of the shot noise matrix £ (see Fig. 6 in [49|). 
We denote these components with a plus- and a minus 
subscript, according to their super- and sub-Poissonian 
shot noise level, respectively. In this manner we can ex- 
pand the clustering signal-to-noise ratio as 



N 

S = b^£-^bP = 



-Po^—P+-^P, (34) 



where Vj and \i are the A^ eigenvectors and eigenvalues 
of the shot noise matrix £ and b± and £± are the effective 
bias and shot noise of the two principal components of 
the tracer density field. 



(35) 



where the summation runs over all individual objects in 
the volume. The w±{M) are weighting functions corre- 
sponding to the two non-Poisson eigenvectors of the shot 
noise matrix, as depicted in Fig. [T] Equation (|34l) states 
that a splitting of the tracer density field into A^ mass 
bins yields, in the limit A^ — > oo, about the same clus- 
tering signal-to-noise ratio as simply considering the two 
weighted fields 5+ and 5^. Therefore, these are the most 
promising candidates to carry out a multitracer analysis 
with. 

In [i^ the functional form of W- (M) was found to be 
well described by 



w_(M) = M + Mo , 



(36) 
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where Mq is about three times the minimum halo mass 
resolved in the simulation, Mq ~ 3M„iin- For the second 
weighting function we find 



w+{M) ^ M ~ Ml + M^/M 



(37) 



to be a good fit to the second eigenvector of the shot noise 
matrix, as shown in Fig. [T] For the constants Mi and M2 
we find Ml ~ 1 X lO^^h'^MQ, M2 ~ 3 x lO^^h'^MQ at 
z = and Ml ~ 1 X IO^/i-^Mq, M2 ~ 1 x IO^^/i-^Mq 
at z = 1. 

The physical origin of the first principal component 
is related to halo exclusion effects [50|. The sampling of 
the density field with halos is less stochastic than Poisson 
sampling (sampling with points) due their finite exten- 
sion. Since the exclusion volume of halos is proportional 
to their mass, this effect is strongest at the high mass end 
(49} . The second principal component can be interpreted 
as a loop correction to the galaxy bias [IH , coming from 
the second-order term in a local bias expansion model 
(52} . Due to its nonlinear character, it adds a super- 
Poissonian shot noise contribution to the two-point clus- 
tering statistics of halos, originating from the squared 
density field [ssl [s^. However, through mode-coupling 
it also yields a second-order clustering signal that orig- 
inates from the bispectrum (three-point function) and 
adds valuable information coming from smaller scales. 

Figure [2] displays the three covariance matrix elements 
C++ = (S+dl), C+- = {S+S*_) and C__ = (<5_(5* ) in 
real-space (/i = 0) at redshift z = 0, extracted from our 
simulation. The three power spectra are obviously highly 
correlated and closely follow the shape of the estimated 
dark matter power spectrum up to fc ~ O.l/iMpc"^. By 
construction, the shot noises of the two fields are not cor- 
related, i.e. (e-e+) = 0. Taking the ratios as in Eq. (fT4| 
yields three possible estimators for the relative galaxy 
bias a = bg_ /bg^ ~ 1.6. In the following two subsections, 
we will provide evidence for the claim that the fields (5+ 
and (5_ are indeed the optimal choice for a multitracer 
analysis. 



2. Sampling variance cancellation 

The idea of utilizing multiple tracers is to cancel sam- 
pling variance from the underlying density field 5. To 
quantify the magnitude of cancellation between any two 
tracers Jg^ and , we define the following statistic: 



2 



(|&2'5gi(k,M) - &i(?g3(k,^)| ) 
(|62<5g,(k,^)|2) + (|Mg.(k,A*)r 



(38) 



Here, hi and 62 is the effective bias as defined in Eq. (jS]). 
If the two tracers are completely uncorrelated (r = 0), 
there is no cancellation and (Tsv — 1- Yet, if they are 
perfectly correlated (r — 1), crgv = 0. Because the real 




0.1 

k [/iMpc 

FIG. 2. Real-space auto- and cross power spectra of the two 
fields S+ and S- (as indicated), obtained through weighting 
the halo catalog by the two non-Poisson eigenvectors of the 
shot noise matrix w+ and W-. The long-dashed black line 
shows the dark matter power spectrum from the simulation 
and the dotted line its linear theory prediction. The three 
different estimators for the relative bias a between S+ and 5- , 
as defined in Eq. (fTl)) . are depicted in dot-dashed {C-\ — /C++), 

dot-dot-dot-dashed (C /C-| ) and dashed (-^C--/C++). 

For visibihty, they were shifted upwards by a factor of 10"^. 



and imaginary part of any given Fourier mode is uncor- 
related, we can swap them for one of the tracers to mimic 
the case of no correlation. 

We consider the following selection criteria for two 
tracers from our halo catalog: 

20/80: lightest 20% vs. heaviest 80% of aU halos. 

50/50: lightest 50% vs. heaviest 50% of ah halos. 

80/20: lightest 80% vs. heaviest 20% of aU halos. 

u/w-i all uniformly weighted vs. all w--weighted halos. 

w+/i(j_: all t(;+-weighted vs. all w_-weighted halos. 

The first three are simply obtained via cutting the halo 
catalog in two at different mass thresholds Afcut to 
yield the indicated abundances in each bin. The fourth 
selection utilizes the whole uniform halo catalog (not 
weighted) and its w_-weighted form, while the last one 
uses both weighting functions w± to construct two trac- 
ers from one and the same halo catalog. Further de- 
tails about the tracers with these selection criteria can 
be found in Tab. ID 

Figure [3] shows ctsv for the different tracer pairs, on 
large scales they all exhibit significant cancellation of 
sampling variance as compared to the reference case with 
switched real and imaginary parts (dotted lines). To- 
wards smaller scales this effect is deteriorated due to the 
onset of nonlinearities and velocity dispersion (FoG ef- 
fects) [2^. As evident from the plot, a combination of 
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TABLE I. Details of the tracer selection 



Selection 




Ml 


M2 


ni 




n2 






£22 


£12 


hi 


bs2 


Smax 












] 




] 






[h-^Mpc^] 








20/80 


1.20 X 10^^ 


1.04 X 10" 


7.72 X 10" 


8.50 X 10" 


4 


3.15 X 10" 


3 


1432 


964 


448 


0.899 


0.988 


74 


50/50 


2.00 X 10" 


1.33 X 10" 


1.13 X 10" 


2.00 X 10" 


3 


2.00 X 10" 


3 


911 


1243 


547 


0.896 


1.042 


73 


80/20 


5.00 X 10" 


1.97 X 10" 


2.22 X 10" 


3.15 X 10" 


3 


8.50 X 10" 


4 


854 


1770 


606 


0.901 


1.217 


73 


u/w- 




6.30 X 10" 
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FIG. 3. Sampling variance statistic asv for 5 different tracer 
selections at z = as described in the text, 20/80 (dot-dot- 
dot-dashed, blue), 50/50 (dot-dashed, green), 80/20 (dashed, 
red), u/w- (solid, orange) and w+jw- (long-dashed, yellow). 
Dotted lines show the corresponding results for uncorrelated 
modes by swapping the real and imaginary parts of one of the 
tracer's Fourier modes. 



the uniform with the -weighted halos yields the high- 
est cancellation of sampling variance. However, the com- 
bination of i(;+-weighted and w_ -weighted halos shows a 
comparable suppression, as opposed to cutting the halo 
catalog in two, which yields less cancellation, especially 
for a low mass cut. If a mass cut is imposed to construct 
two tracers, the highest sampling variance cancellation 
is achieved when the same abundance of objects in each 
of the resulting catalogs is chosen (50/50), which is in 
agreement with the findings of [ssj . 



3. Fisher information 

Cancellation of sampling variance alone is not a suffi- 
cient indicator on how well cosmological parameters can 
be constrained in a multitracer analysis. This is because 
one is looking for relative changes in the clustering signal 
from multiple tracers, and not for the absolute clustering 
amplitude in each tracer. In this paragraph we will show 
that it is desirable to have tracers with a high relative 



galaxy bias ratio a. 

Quantitatively, the achievable accuracy on a given cos- 
mological parameter 9 is determined by the inversion of 
the Fisher matrix of Eq. (HI]). For the sake of simplicity, 
let us consider only the /-/-element. 



s2 + s (ee^^ 



y2 



l-hE2 



(39) 



In the high signal-to-noise regime, E ^ 1, and express- 
ing E, E/, as well as E// in terms of the two principal 
components as in Eq. (|34|) . we get 



(b--b+) 



(40) 

As a second approximation, we can make use of the fact 
that b'i/S^ > bl/£+ and f_ < £+ (see Tab. IJ), which 
yields 



1 - 



1 - 



l + /?Ai- 



2 N 2 



'£7 



(41) 



The single-tracer term (as derived in j]36]) is dominated 
by the first principal component of E, while the multi- 
tracer term depends on the bias ratio a = bg_ /bg^ of 
both principal components and the shot noise f + of the 
second principle component. In order to maximize P//, 
it is thus desirable to have a large a and a low at the 
same time. In the special case of uniform Poisson shot 
noise, Eq. ([M)) reproduces the expression derived in j57i] . 

If the dark matter density field is known separately, 
this Fisher matrix element becomes 



P: 



/i^P 



(42) 



thus independent of the effective bias and limited only by 
the low shot noise level of the first principal component. 

Figure m depicts Fff{k,^ = 1) for all of our five tracer 
selections. For each case we further distinguish between 
the following scenarios: 



8 



20/80 


50/50 


80/20 




u/w_ 


— ^ w+/w_ 

\ 




y' " \ 


\ 

\ 


— ■ " 


\ 

. \ 


\ 

^ " ~ ^ 




.^^~~^~^^^\.\., 














. . • 









0.01 0.1 1 0.01 0.1 1 0.01 0.1 1 0.01 0.1 1 0.01 0.1 1 

k [/jMpc '] k [/iMpc '] k [/jMpc ^] k [/iMpc '1 k [/iMpc '] 



FIG. 4. Fisher matrix element F// as a function of fc at ^ = 1 and 2 = for various tracer selections (indicated in the top 
right of each panel). Results are shown for a single-tracer analysis with the first (solid, blue) and the second tracer (dotted, 
green) used separately, both of the tracers combined in a multitracer analysis (dashed, red) , and each of the tracers combined 
with idealistic dark matter observations (dot-dashed orange and long-dashed yellow, respectively). 



• single-tracer analysis with each of the two tracers 
taken separately 

• multiple-tracer analysis with both tracers combined 

• combined analysis of each tracer with the dark mat- 
ter density field 

Clearly, in a single-tracer analysis the tracer with the 
lowest galaxy bias yields the highest information on the 
growth rate /. This is evident from Eq. p3|. where 6g 
enters in the denominator of the first term and the shot 
noise £ is negligible if £/h'^P <^ 1. Since bg^ is very 
similar in all five cases, the Fisher information from a 
single tracer cannot be increased much by any particular 
choice of tracer. 

On the contrary, a combination of two tracers can can- 
cel out sampling variance and therefore considerably in- 
crease the available information on / if the tracers are 
selected appropriately. As evident from the dashed red 
lines in Fig. 01 the multitracer term in Eq. (HTI) gains im- 
portance over the single-tracer term as moving from the 
left to the right panel. Again, the two tracers obtained by 
a low mass cut yield the worst results, owing to the fact 
that the shot noise of both tracers is super- Poissonian in 
this selection. It is more optimal to impose a high mass 
cut in order to benefit from the sub- Poissonian shot noise 
level of the heaviest halos [49]. However, the highest 
Fisher information content is obtained when correlating 
the two orthogonally weighted fields 6+ and 5-. The 
main reason for this is the large relative galaxy bias be- 
tween the two fields and their relatively low shot noise 
level (see Tab. Hand Eg. (|IT|)). 

We have also explored the possibility of splitting the 
halo catalog into more than two mass bins, considering up 
to TV = 10 tracers. We find that the Fisher information 
from multiple tracers increases with the number of bins N 
and approaches the result obtained with the two weighted 
fields S+ and 6- in the limit of high N. 



Finally, adding the information from the dark matter 
density field to each one of the tracers increases the in- 
formation content on /. In this case, sampling variance 
inherent in the density field S is known and can thus 
be removed from the halo fields directly. According to 
Eq. (|42l) , the Fisher information is inversely proportional 
to the shot noise of the tracer, so the lowest stochasticity 
weight W- yields the best results. 

C. Multitracer fit 

So far we have investigated the Fisher information con- 
tent on the growth rate using multiple tracers of the LSS. 
The question of how to actually constrain parameters of 
interest from a dataset in the optimal way will be an- 
swered in this section. For this task we want to maximize 
the likelihood function from Eq. (|17p , which is equivalent 
to minimizing its negative logarithm, the chi-square 

X'^J2 A*)C^Mg(k, /i) + ^ In (det C) . (43) 

k 

Here, the covariance matrix C is given by the clustering 
model of Eq. ([7]) and we sum over all individual Fourier 
modes <5g(k, /i) from our halo catalog. When we add the 
dark matter density field as an observable, we use the 
model from Eq. ((T5|) and 6 = {S, Sg) as our data vector. 

Figure [S] presents the fitting results for the RSD pa- 
rameter /3 and the product of growth rate / with the 
normalization of the power spectrum erg from our halo 
catalogs at redshift 2 = 0. While (3 has been obtained 
from a single-parameter fit, we have marginalized over 
the galaxy bias of each tracer as a free parameter in the 
fit for /(Tg . In the left column, the standard single-tracer 
analysis utilizing all objects in the halo catalog is per- 
formed. The best fits along with their Icr-error bars 
are shown both in /c-bins (points with error bars), as 
well as cumulative as a function of A; = /cmax with fixed 




FIG. 5. Fit for the redshift-space distortion parameter /3 (top) and the product fas (bottom) from the two-point clustering 
statistics of halos in an A^-body simulation with effective volume Kff — 6.6h~^Gpc^ and halo-mass resolution Mmin — 9.4 x 
lO^^/i~^M0 at z = 0. LEFT; Conventional single-tracer analysis utilizing all halos (not weighted) from the same catalog. 
RIGHT: Multitracer analysis with the two fields 5+ and obtained through weighting the halo catalog with its principal 
components u;+ and w_ . The fitting results are shown in logarithmic bins of k (points with 1-a error bars) , as well as cumulative 
as a function of fc = fcmax with fixed femin = 0.0048/iMpc~'^ (red solid lines with shaded region). Dotted lines show the linear 
theory prediction with / = f2m^^ and ag — 0.81, dashed lines the cumulative sampling variance limit. 



fcmin = 0.0048ft.Mpc~^ (solid lines with shaded region). 
The constraints on /? and fas are clearly affected by sam- 
pling variance, as evident from the large scatter of the 
points at low k. When the number of available Fourier 
modes grows towards higher fc, this scatter becomes 
smaller, however, beyond a scale of fc ~ 0.2/iMpc~^, lin- 
ear theory breaks down and the fits depart from their 
scale-independent linear value assuming / — il^j^^ (dot- 
ted line) . The cumulative sampling variance limit for the 
determination of fag is shown as a dashed line. Clearly, 
the single-tracer fit yields a substantially larger uncer- 
tainty compared to this limit. 

When combining the two tx;+/w_ -weighted tracers in 
a multitracer analysis as shown in the right column of 
Fig- El the scatter of the fit is significantly suppressed. On 
the largest scales, the errors are reduced by up to a factor 
of 4 and the constraints on /erg reach the sampling vari- 
ance limit closely. However, the fit seems to deviate from 
the linear theory prediction already at fc ~ 0.04/iMpc~^. 
On the other hand, the apparently more linear behavior 
of the single-tracer analysis may likely be coincidental at 
k > 0.04/iMpc~^. This is supported by results shown in 



the left panel of Fig. [51 Here, we combine the uniform 
halo catalog (without weighting) with the dark matter 
field from our simulation to fit for / directly. Obviously, 
sampling variance has decreased even further, but devi- 
ations from linear theory already kick in at a scale of 
k ~ 0.04/iMpc"\ 

Most impressive constraints on the growth rate are ob- 
tained when we combine the w_ -weighted halos with the 
dark matter density field, as depicted in the right panel 
of Fig. [6l In this case, sampling variance has almost 
canceled out completely and the error bars on / have 
diminished by up to a factor of 142 (at the peak of the 
power spectrum) , when compared to the standard single- 
tracer analysis. Moreover, deviations from linear theory 
are very small up to scales of fc ^ 0.2ft,Mpc~^, making 
this kind of experiment the most promising one out of 
the four considered scenarios. 

Methods to combine galaxy clustering and weak lens- 
ing data have been studied extensively in the recent 
literature (e.g., [58l - [63 |). suggesting high improvements 
for precision cosmology. Unfortunately, obtaining the 
dark matter density field in 3D from observations is a 




FIG. 6. Fit for the growth rate / from the two-point clustering statistics of halos and the dark matter combined. LEFT: All 
uniform halos (not weighted) and the dark matter. RIGHT: All halos, weighted with the lowest stochasticity weight w^, and 
the dark matter. The meaning of lines and symbols is the same as in Fig. [5] 



highly nontrivial problem and is subject of active re- 
search. Weak lensing tomography is the technique aiming 
to achieve this goal (e.g., ,65n6&J), but the resolution in 
the radial direction is not expected to be high, because 
of the relatively broad lensing kernels along the line of 
sight. Moreover, we assume an ideal reconstruction of 
the dark matter density field without considering addi- 
tional sources of error involved in the lensing measure- 
ment, such as shape noise and intrinsic alignment^ for 
example. 

Without knowledge about (5, in principle we would have 
to marginalize over all the parameters O'^^^ of the dark 
matter power spectrum as well. This implies calculating 
the transfer function in each iteration of the fitting proce- 
dure, which goes beyond the scope of this paper. In the 
high signal-to-noise regime of the multitracer analysis, 
the degeneracy between the parameters O'^^^ and O'^^^ is 
expected to be rather weak (except the fundamental de- 
generacy between / and cts), as the mixed terms between 
E,; and Pj in Eq. (1^ are suppressed by S. Moreover, in 
Appendix C of [51] it has been shown that P{k) cancels 
out to a high degree in the chi-square of Eq. P5|) . Of 
course, in case the dark matter density field is available, 
we do not have to worry about those issues, since 6g and 
P are directly observable. 

In order to quantify the gains in accuracy, we com- 
pare the size of the error bars on j3 and /cts from the 
multitracer- to the single-tracer analysis in Fig. [7] (solid 
and dashed blue lines, respectively). Obviously, sampling 
variance mostly cancels on the largest scales, yielding 
improvements in accuracy of up to a factor of 4. Be- 
yond scales of fc ~ 0.1/iMpc~^, the improvement is de- 
teriorated due to the onset of nonlinear clustering and 
mode-coupling. Therefore, the drop of the curves at 
k > 0.3ft.Mpc~^ is likely an artifact and should not be 
trusted. At redshift z = 1 this turnover is moved to 



smaller scales, the overall improvement in the error ratio 
is however slightly deteriorated. In order to access the 
cosmological information content buried in the semilinear 
regime of galaxy clustering in redshift-space, one can not 
avoid having to invoke more elaborate models involving 
perturbative methods, such as the ones proposed in [70| . 

More gains can be achieved when galaxies and ideal 
dark matter observations are combined, the improve- 
ment in the accuracy on j3 compared to the ordinary 
single-tracer analysis amounts to about a factor of 10 
in this case (solid green line). If, additionally, galaxies 
are weighted optimally, it increases by another factor of 
10, two orders of magnitude better than what a single- 
tracer analysis can achieve. In this case the improvement 
even persists down to smaller scales of fc ~ 0.3/iMpc^^. 
However, the effect of optimal weighting is diminished 
towards higher redshifts. 

Unfortunately, the halo masses we used to construct 
our weighted density fields are not directly observable in 
reality. Yet, they correlate with many observables, such 
as X-ray luminosity, galaxy richness, weak lensing shear, 
velocity dispersion or the thermal Sunyaev-Zel'dovich 
(SZ) effect (e.g., (Zli)- Scaling relations between these 
observables and halo mass can be calibrated with nu- 
merical simulations to obtain unbiased mass-proxies with 
minimal scatter (e.g., [zSjIzJ)- While optical methods 
show a rather large scatter of uin a/ — 0.45 [T^I, X-ray or 
SZ observations yield tighter relations with uinM reach- 
ing below 0.1 [zsj. 

We artificially add a constant log-normal scatter to 
the halo masses of our catalog in order to mimic the 
observational uncertainties in the mass determination 
for a rather pessimistic {uxaM = 0.5) and a more op- 
timistic scenario {oinM = 0.1). The results are depicted 
in Fig. [7] as dotted and dot-dashed lines, respectively. 
While the constraints from the multitracer analysis on 
/3 are only marginally affected for (T\^m = 0.1 and de- 



11 




FIG. 7. Ratios of the binned one-sigma error bars (from Fig. [5]) on /3 (solid blue) and fas (dashed blue) when comparing the 
optimal multitracer- to the single-tracer analysis at 2; = (left panel) and z = 1 (right panel). Additionally, the improvement 
in constraints on /3 from a combined clustering analysis of halos and dark matter as compared to the single-tracer case is shown 
in solid green (no weighting of halos) and in solid red (optimally weighted halos). Adding a log-normal scatter of (JinM = 0.1 
(dot-dashed) and crinM = 0.5 (dotted) to the halo masses results in a degradation of the constraints from the weighted fields. 



grade by roughly 20-30 percent for (Tin a/ = 0.5, the com- 
bined analysis using optimally weighted halos and the 
dark matter is more severely deteriorated by mass scat- 
ter. Here, even the optimistic scenario increases the un- 
certainty on the growth rate by roughly 50 percent, while 
the benefits from optimal weighting are completely lost 
when going to trinM = 0.5. Clearly, the high level of pre- 
cision that can be obtained with this kind of experiment 
demands precise mass estimates. Fortunately, a whole 
industry of existing and planned experiments devoted to 
cluster cosmology will provide those high quality data 
(e.g., [Bill). 



D. Halo model predictions 

In this section we want to investigate how our results 
depend on the resolution of the simulation, respectively 
the minimum resolved halo mass. In a real experiment, 
this corresponds to the depth of the galaxy survey, with 
a corresponding luminosity threshold. Because iV-body 
simulations with a given volume become increasingly ex- 
pensive with higher resolution, we will turn to theoretical 
predictions henceforth. 

The clustering properties of dark matter and halos can 
be neatly described by the halo model (see, e.g., [s^). 
The basic idea is to separately describe the clustering 
within a given halo (one-halo term) and the clustering 
amongst different halos (two-halo term) . In [4^ the halo 
model is utilized to derive an analytical expression for 



the shot noise matrix. The result can be written as 

e = n-^I - 6gA4t - Mbl , (44) 

where n is the number density of halos per bin, I the 
identity matrix, M = M/pm - b^{nM'^)/2p^ and M 
a vector containing the mean halo mass of each bin. In 
appendix [B] we utilize this expression to derive the clus- 
tering signal-to- noise ratio E (as well as and ) from 
the halo model. With Eq. (15^ we can then determine 
the Fisher information on the growth rate and compare 
the single-tracer to the multitracer analysis. 

Figure |5] displays the same ratio as Fig. [7] for the un- 
certainty on P, but now as a function of minimum halo 
mass M„ii,i at a fixed scale of A; ~ 0.016/iMpc~^ (peak 
of the power spectrum). The gains from the multitracer 
method kick in at Afniin — lO^'^/i^^M©, where S ~ 1, and 
increase towards lower Af„iin due to the growing signal- 
to-noise ratio. A combination of galaxy- and dark matter 
observations may increase these gains further, especially 
when the galaxies are weighted optimally. In this case 
there is no saturation towards lower Mmin and the im- 
provement compared to the single-tracer analysis con- 
tinues to grow. In contrast, the uniform galaxy over- 
density field (not weighted) combined with dark matter 
already shows up a saturation at Mmin — 10^°/i~^Mq, so 
no more information on the growth rate can be gained 
when even more lighter halos are included in this kind 
of analysis. Also note that the halo model predictions 
underestimate the improvements obtained when adding 
dark matter clustering information, as our A^-body re- 
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suits with Minin ^ 9.4 X 10"/i"^Mq yield higher im- 
provements (see Fig. [7]). 



IV. CONCLUSIONS 

In this paper we investigated the benefits of using 
weights in a multitracer analysis of LSS, with a partic- 
ular focus on constraining the growth rate of structure 
formation. On the basis of earlier results on the cluster- 
ing properties of dark matter halos and their stochasticity 
|49| , we argue that the gains from a multitracer analysis 
in the sense of [l^l can be achieved by considering only 
the two principal components of the clustering signal- 
to-noise ratio E (or, equivalently, the two non-Poisson 
eigenvectors of the shot noise matrix £). 

We present their explicit functional forms in terms of 
weights, showing that the first one coincides with the 
weighting function explored in previous work j49i] , giving 
rise to low stochasticity and high bias. For the second 
one the weights are also mass dependent, but have a zero 
crossing, such that the overall bias is low. This yields 
a high relative galaxy bias a between the two tracers, 
maximizing the Fisher information content on the cos- 
mological parameters [l3] ■ All of the other eigenvectors 
oscillate around zero and add very little information. The 
advantage of reducing the information to two eigenvec- 
tors is that all of the objects in a given catalog can be 
used to construct the two principal components, while 
in the conventional multitracer analysis the catalog has 
to be split into many lower number density sub-samples 
with higher shot noise [s^ [s^] . 

On the basis of numerical TV-body simulations of dark 
matter halos, we demonstrate that the constraints on /3 
and /(Ts can be improved by up to a factor of 4 rela- 
tive to a single-tracer method, but most of the improve- 
ment comes from large scales (low fc), while for higher k 
the gains are smaller and vanish above k ~ 0.1/iMpc^^, 
where nonlinear effects introduce additional stochasticity 
between the two tracers. This technique is fairly insensi- 
tive to the observational uncertainty on the halo masses, 
as even a 50% log-normal scatter does not degrade the 
improvements significantly. Halo model considerations 
suggest even higher gains of the method with increasing 
mass resolution. 

One potential concern for our method is the possibility 
that galaxies might be bad tracers of their host-halo cen- 
ters |8j] and therefore exhibit less pronounced principal 
components in the clustering signal-to-noise ratio that 
are distinct from Poisson sampling. However, there are 
strong indications that certain types of galaxies do show 
strong correlations in both position and mass with their 
host halo, e.g. luminous red galaxies (LRGs) [s^. Tech- 
niques to distinguish satellite- from central galaxies have 
been developed and the satellite-fraction can be used as 
an estimator of the host- halo mass [s^ . Mock LRG cata- 
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FIG. 8. Halo model prediction for the improvement on 
from the single-tracer to the optimal multitracer analysis 
(blue) at fc ~ 0.016/iMpc~^ as a function of the lower halo- 
mass threshold Mmin at z = (solid) and z = 1 (dashed). 
The results for the combined analysis of uniform halos (not 
weighted) and dark matter (green), as well as the combina- 
tion of optimally weighted halos with the dark matter (red) 
are also depicted. 



logs obtained from a halo occupation distribution (HOD) 
model suggest some reduction in stochasticity is possible 
even without explicit knowledge of halo masses [8g| . 

Whether these gains translate into a useful constraint 
on the final cosmological parameters depends on our abil- 
ity to model nonlinear RSD effects. We find nonlinear 
effects are important for /3 already at fc > 0.03/iMpc~^, 
although they appear to be important for /ag only at 
k > 0.1/iMpc~^. In the most pessimistic case where 
the RSD model cannot be trusted for k > 0.03/iMpc~^, 
the multitracer method provides major gains relative to 
the single-tracer case, but neither method provides very 
strong constraints overall, because of the limited num- 
ber of available Fourier modes. In the case where we 
can use all the modes up to fc ~ 0.1/iMpc~^, the over- 
all errors are considerably smaller and the multitracer 
method provides less of an advantage. It is clear that a 
better modeling of the nonlinear effects in RSD is needed 
to understand the ultimate reach of RSD in both single- 
and multitracer methods. 

In a more idealistic scenario, we also consider the joint 
analysis of halos and the dark matter density field, which 
in principle is achievable via a combination of spectro- 
scopic redshift surveys and weak lensing tomography. 
Here, utilizing optimal weights can yield up to 2 orders of 
magnitude improvements in constraining /3 as compared 
to a single-tracer analysis, but the method is more prone 
to uncertainties in the halo mass estimates. It is unlikely 
that this gain can be achieved in practice, since it is very 
difficult to measure dark matter clustering in the radial 
direction directly. 
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A further technique to construct differently biased 
tracers of the density field makes use of nonhnear trans- 
formations [S^]- Although it is difficuh to describe the 
effects of a nonhnear transformation on both signal and 
noise in galaxy clustering data, combined with optimal 
weights this may provide another tool for the multitracer 
analysis. 

In this paper we have focused on the information that 
can be extracted from RSD, in particular /3 and /cts, but 
our method is not limited to constraints on the growth 
rate, but may be applied to the analysis of primordial 
non-Gaussianity [5l[, general relativistic corrections in 
large-scale clustering [s^] , the Alcock-Paczyhski test [l2| , 
or any other quantity that influences the effective bias of 
tracers of the density field. It is possible that a better 
model of nonlinear RSD may yield a more efficient mul- 
titracer method, where the gains relative to the single- 
tracer analysis described here on large scales can be ex- 



tended to smaller scales. We leave these directions for 
the future. 
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Appendix A: FISHER MATRIX FOR MULTIPLE BIASED TRACERS 



With Eqs. (Uni) and ^ plugged into Eq. (HH]), we have 
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With additional knowledge about the dark matter density field we need to work out Eq. by plugging in Eqs. (pS)) 
and (P71) . Let us first note that 
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Appendix B: HALO MODEL PREDICTION FOR THE CLUSTERING SIGNAL-TO-NOISE RATIO 

In the halo model the shot noise matrix is given by Eq. (|44p . In order to invert 5, we write £ = A. A4bg with 
A = fi^^I — bgA^^ and apply the Sherman-Morrison formula twice: 

= + ^— , A-^ = nl + § ^ . Bl 

1-blA-^M n-i-M^bg 

With 

blA-^b, (l - blA-^Ai) + blA-^MblA-^b, 

E,- ^ bl£-%P = ° I ' , (B2) 

1 - blA M 



and 



bib, (n-i - M^b,) + bjb.M^b, 
blAr% = ^ . , n , (B3) 

1 - blA-'M = A y ' n , (B4) 



and similar terms combining b, bi, bg and A^, after some algebra we get 



bjMblbgM^bj + (blbgM''b, + bjAib^bJ - M^'bA + blb^M^Mblb, 
= blbjfiP + ^ ^^-^ nP . (B5) 

In the continuous limit (N — > oo), we can exchange the vector products by integrals over the mass function and set 
^tot = nN. This finally yields 

^ ^ {bl){Mh){Mb,) + ((b.bg)(>(b,) + {Mh){b,bg)) (ntol - (.Mbg)) + {b\){b,bg){M') _ 

{n^:,~{Mbg)f-{bl){M^) 

(B6) 
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